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Abstract 



We present predictions for the radio pulses emitted by extensive air showers us- 
"^ I ing ZHAireS, an AIRES-based Monte Carlo code that takes into account the full 

complexity of ultra-high energy cosmic-ray induced shower development in the at- 
mosphere, and allows the calculation of the electric field in both the time and 
"^ I frequency domains. We do not presuppose any emission mechanism and our re- 

^; suits are compatible with a superposition of geomagnetic and charge excess radio 

emission effects. We investigate the polarization of the electric field as well as the 
^N I effects of the refractive index n and shower geometry on the radio pulses. We show 

Qv ■ that geometry, coupled to the relativistic effects that appear when using a realistic 

OO ■ refractive index n > 1, play a prominent role on the radio emission of air showers. 

r-^ ■ Key words: high energy cosmic rays and neutrinos, high energy showers, 

^D I Cherenkov radio emission 

PACS: 95.85.Bh, 95.85.Ry, 29.40.-n, 



C^ ■ 1 Introduction 



In the last decades, the study of ultra high energy cosmic rays (UHECR) has 
been one of the most active areas in astroparticle physics [Tj. The detection 
of extensive air showers (EAS) created by UHECR in the atmosphere has 
been accomplished mainly with two detection methods. The first consists on 
detecting the particles of the cascade reaching the ground using an array of 
particle detectors. The second method uses telescopes to detect the fluores- 
cence photons emitted by the particles in the shower as they travel through 
the atmosphere. Only the latter method is able to directly measure the EAS 
longitudinal development in the atmosphere, but is subject to a very low duty 
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cycle (~ 10%), since it can only be used in clear moonless nights. The surface 
array and fluorescence techniques are used simultaneously at the Pierre Auger 
cosmic ray observatory |2]. 

Radio emission from EAS was first observed by Jelley et al. |[3j in 1964. The- 
oretical and experimental research in this area was very active in the sixties 
and early seventies. We refer the reader to a review by Allan |1] and references 
therein. Nevertheless, interest in this technique declined in the late seventies, 
mainly due to radio interference [5]. Recent developments in high speed elec- 
tronics and information technology have renewed interest in the detection of 
radio emission in the MHz range from the particles in EAS. The radio de- 
tection technique is, in principle, sensitive to the longitudinal development of 
the shower, like the Fluorescence technique, but its duty cycle is much higher, 
since it could be used anytime except during thunderstorms. 

This resurgence of the radio technique has taken form as new experiments 
have been developed, such as CODALEMA [6], LOPES [5)17118] and AERA |9], 
accompanied by new calculations of the radio emission in EAS, which in- 
clude analytical techniques with different levels of sophistication 
Monte Carlo methods J14|15lll6lll7fl8fl9] and semi-analytical methods 

In this work we present ZHAireS [2TIT9] (ZHS-I- AIRES), a new simulation of 
radio emission in EAS that combines the full shower simulation capabilities of 
AIRES |22j with specific algorithms developed to calculate the electric field 
emitted by particles in dense media showers, implemented in the well tested 
ZHS code [23|2l] . These algorithms are obtained from first principles, so no 
emission mechanism or model is presupposed. The radio emission calculations 
are done in parallel to the AIRES shower simulation: As each charged particle 
in the shower is propagated by AIRES in steps, each propagation step is taken 
as a single particle track and its contribution to the radio emission is calculated 
and added to the total electric field in both the time and frequency domains. 
This procedure naturally accounts for interference effects associated to the 
different space-time positions of the particles in the shower. 

Other Monte Carlo simulations of radio emission of air showers exist, such 
as REAS [TBPT] , where the electron and positron tracks from CORSIKA [SS] 
simulations are first histogrammed and then used to generate random e^ tra- 
jectories, which in turn are used for the radio emission calculations in the time 
domain only. An earlier version of this code, REAS2 [16], was based only on 
geosynchrotron emission, but was shown to be inconsistent with later simula- 
tions [TTfT^ . which in turn are based on a pretty generic algorithm [3] that has 
been used for a long time to simulate pulses generated in dense media [23|2^ . 
Another difference is that ZHAireS uses a model for the variation of the re- 
fractive index with altitude, while REAS uses a fixed refractive index equal 
to unity in its calculations. 



This paper is organized as follows: In section |2] we give a short overview of the 
main radio emission mechanisms in air showers and the specific formalism used 
for emission calculations in ZHAireS, including a model for the variation of the 
refractive index with altitude. In section [3] we develop a one- dimensional toy 
model useful for understanding the main characteristics of the emission and 
stress the prominent role played by geometry. In section H] we show results of 
ZHAireS simulations, discussing the influence of the distance from the antenna 
to the shower core, the refractive index and the shower zenith angle on the 
electric field pulse. We also discuss the spectrum of the radio emission. In 
section \5\ we analyze the polarization of the electric field and discuss how it 
can be used to separate the contributions from various emission mechanisms. 
Finally, in section |6] we conclude the paper. 



2 Radio Emission in Air Showers 



2.1 Main emission mechanisms 



In air showers, the dominant mechanism responsible for the radio emission is 
believed to be the deflection by the geomagnetic field of electrons and positrons 
in the shower [261171127] . Several approaches have been developed to model 
this emission. Since the direction of the Lorentz force depends on charge, it 
leads to a spatial separation of electrons and positrons in the shower, that 
can be thought of as a moving macroscopic dipole and a transverse current 
traveling through the atmosphere at a speed v ^ c along with the shower 
front, and which has been the basis for macroscopic calculations [28|llf20] . 
Another approach, which is adopted in this work, is to calculate the emission 
for each particle trajectory, i.e. a microscopic approach P5|24f21fl9fl7] . 

The geomagnetic mechanism has a quite clear signature. Its polarization is 
anti-parallel to the direction of the Lorentz force, i.e. in the direction of — /3 x 
B, where /3 is the speed of the particle in c units and B the geomagnetic 
field I2811291I301I3T]. 



Another emission mechanism also thought to be important in EAS is the 
Askaryan effect. It was first proposed by Askaryan [3lj, who suggested that 
coherent radiation could be emitted by showers in which a charge excess devel- 
ops. This mechanism, which dominates the radio emission of showers in dense 
media, is also known as the charge excess mechanism. For a particle following 
a straight track at constant speed, it can be shown [32] that: 



Erad oc -e[u X (m X (3)] (1) 



ti,Ei,x 




(const) 



Fig. 1. Scheme of a single particle track. Speed and direction to the observer are 
assumed constant inside each track. 

So the electric field emitted by particles traveling along the shower axis lies in 
the plane defined by the direction of the shower axis (/3) and the observation 
direction u. Furthermore, it is perpendicular to u and its direction depends 
on the charge of the particle, pointing towards (away from) the shower axis 
for negative (positive) charges [23|2^ . This means that if there is no charge 
excess in the shower, the net electric field is zero. However it is well known 
that knock-on interactions and Compton scattering incorporate electrons from 
molecules of the medium into the shower, leading to an excess of moving 
negative charges [34J , which in turn are responsible for a net electric field with 
a radial polarization w.r.t. the shower core. 



2.2 Radio emission calculation in the ZHAire S code 



The algorithms used for the calculation of radio emission in ZHAireS are based 
on ZHS algorithms [23p^ . These were derived from first principles, namely 
from the Lienard-Wiechert potentials, and thus do not assume any emission 
mechanism. The full derivation can be seen in [23j for frequency domain cal- 
culations, and [23] for the time domain. The trajectories of shower particles 
are divided in tracks, which can be made arbitrarily small. Both the speed of 
the particle and the direction to the observer are assumed not to vary over 
the track. In Fig. [1] we show a schematic picture of such a track. A conve- 
nient division of trajectories in small tracks is already performed by Monte 
Carlo simulations of shower development in order to propagate the particles 
in the shower. The ZHS algorithm is then used to calculate the contribution of 
each track to the net emission of the shower, accounting for any interference 
between tracks. This approach automatically takes into account the contribu- 
tion to the electric field due to the start, end, and any change in the direction 
and energy of each particle track. Thus any kind of deflection, scattering, cre- 
ation or annihilation of a charged particle, due to any physical process used 
to simulate the shower is taken into account in the radio emissionlT9]|21ll37l. 



In the original ZHS algorithms [23)12^ . the field calculation is performed in 
the Fraunhofer approximation. In this case the observer "sees" the radiation 



arriving from all the points in the shower at the same angle, and only changes 
in the phase of the contribution from each point are taken into account through 
a simple projection. This approach works well for showers in dense media, such 
as ice, since the size of the shower is much smaller than the typical distance 
to the antennas. But in the case of air showers, the distance to the observer is 
usually of the order of the shower size, and thus the Fraunhofer approximation 
breaks down. In order to make the method valid for the closer observers in 
air showers, we allow the distance and direction to the observer to change 
from one track to the next. The direction and distance to the antenna are 
still taken to be constant inside a single track (see Fig. [T]), i.e. we still use 
the Fraunhofer approximation inside each track. It can be shown [36] that 
this procedure reproduces the expected behavior E oc l/y/R of the field [38] . 
where R is the distance between the charge and the observer. Since we still use 
the Fraunhofer approximation inside a single track of length L, the following 
condition has to be satisfied: 



L^ sin^ e A 

< — (2) 

where A is the wavelength of the emission. In the simulation, the vast majority 
of tracks satisfy this condition for frequencies up to 300 MHz, which is much 
higher than the frequency of the maximum of the emission, ~ 1 — 30 MHz. 
If a particular track in the simulation does not satisfy this condition, e.g. a 
track passes very close to an antenna, it is further divided into several sub- 
tracks, and the field calculation is performed using a different R and 6 for each 
sub-track. 

The positions Zi 2, times ti 2 and kinetic energies Ei^2 for the beginning and 
end points of the track are obtained directly from AIRES. This, along with 
the position Xant of the antenna is all that is needed for the calculation of the 
contribution to the vector potential A{t,u) due to a given track: 

AtiRc l-n/3 -u 

where u = Ru is the vector from the middle point of the track to the antenna, 
/3 = -u/c, /3^ = — M X (n X /3) is the projection of (5 onto a plane perpendic- 
ular to M, tf^l = ^1,2 + nR/c — nf3 ■ u (ti_2 — ^0) are the retarded (detection) 
times for the beginning and end of the track, respectively, to = (^1 +_^)/2 is 
the average time for the track and 6(x) is the Heaviside step functioriU. Note 
that Eq. (^ is written in the radiation gauge, since we disregard the static 
term of the Lienard-Wiechert potentials in its derivation 



^ In the formalism, as described in ^24j, the limit of eq. ([3]) for (1 — n/3 • -u) — t- is 
used instead of eq. ([3]) if the denominator vanishes. 



To obtain the radio signal at each antenna for the shower as a whole, the 
contributions of each particle track to the vector potential A{t) are added 
upLil. This net vector potential is then differentiated with respect to time to 
obtain the net electric field as a function of time for each antenna. In the far 
field, this formalism is equivalent to the one used in [17], as shown in 



Besides the field calculations in the time domain, ZHAireS can calculate the 
radio emission in the frequency domain as well. The radiation term of the elec- 
tric field is E(t) = —dA/dt. Applying this to the vector potential expression 
(eq. [3]) we obtain 



g(,,) = -^;?/('-'f')-5-'"') (4) 



In the ZHS formalism [23|2l] . we use the following convention for the Fourier 
transform: 

_ poo 

fico) = 2 / f{t)e'-'dt (5) 

Applying this Fourier transform convention to Eq. (j4]) we obtain [23] : 

E{uj, u) = ^^pj ^ (6) 

which is the expression we use for the Fresnel regime frequency-domain cal- 
culations in ZHAireS. The contributions of each track to E{u) are added up 
to obtain the net spectrum at each antenna. 



2.3 Variable refractive index 



As stated before, in ZHAireS we use a model for the variation of the refractive 
index n with altitude in the atmosphere. In fact there is a strong dependence 
of n on temperature, pressure and humidity. To take the dependence with 
altitude into account we use an exponential model for the variation of the 
refractivity 7^ with height h, motivated by the exponential decrease of density 
with altitude, 

n{h)=nsexp{-Krh) (7) 

where TZ{h) = [n{h) — 1] x 10^ is the refractivity at an altitude h in km and we 
used TZs = Tiih = 0) = 325 and Ky = 0.1218 km~^. These values reproduce 
the refractivity calculated in [39] up to /i ~ 10 km within less than ~ 1%. The 
values in [39] take into account the humidity dependence of the refractivity. 



^ Note that interference effects are taken into account automatically, since the con- 
tribution of each track to the vector potential is related to a specific retarded time. 



which increases the refractivity at low altitudes where the effect of a variable 
index of refraction is most important, since the shower has a larger number of 
particles. At higher altitudes, the exponential model slightly overestimates the 
refractivity, but the effect of the variable index of refraction at these altitudes 
is much less important, since above 20 km the shower has barely started 
developing and the number of particles is small. It is also worth noting that 
since TZ depends on temperature and humidity, there are large seasonal and 
even daily variations in the refractivity that are associated to atmospheric 
conditions and cannot thus be accurately described by any modetil. 

The variable atmospheric refractive index would in principle make the radio 
emission follow a curved path. However, in [20] it was shown that the effect 
of the deviation from a straight path on the time structure of the pulse is 
negligible, and thus in ZHAireS we assume that the radio emission follows a 
straight path from the emission point to the antenna. However, we explicitely 
take into account the effect of a variable refractive index in the propagation 
time of the signal. For this purpose we calculate an effective refractive index 
rie// for each particle track in the ZHAireS shower simulation: 



nejf = 1 + 7^e// X lo-^ 7^e// = -^ j^ n{h) di (8) 



where R is the distance from the track to the observer, and dl is an infinitesimal 
length along that path whose altitude h varies along the path. It is important 
to note that while rie// is used for the calculation of the retarded times, n{h) 
is used for the angular dependence of the emission, given by the denominator 
of Eqs. ([3]) and (E]), i.e. the Cherenkov angle depends only on the refractive 
index at the emission altitude. 

We expect the effect of the variable refractive index to be more important at 
relatively low altitudes, where n is larger. Since n in our model varies between 
n ~ 1 at high altitudes and n ~ 1.000325 close to ground, we expect the 
effect of the variable n on the pulse characteristics discussed in section [3] (such 
as start-time, duration and compression in time) to be between the effects 
obtained for n = 1 and n = 1.0003. This can be clearly seen in Figs. OandE] 
below. 



^ If the humidity as a function of height can be monitored using e.g. a LIDAR, then 
more accurate values of the refractivity for specific atmospheric conditions could be 
calculated by adding a humidity term to Eq. ([Tj) , as described in |39j . 




4. '_A J- antenna 

Fig. 2. Sketch of the one-dimensional model. The shower development is represented 
by a ID line neglecting the lateral structure. The shower front is assumed to travel 
at a speed c and reach ground at a time t' = 0. Radio emission from a height h 
travels along a distance R before arriving at the observer located at a distance r to 
the shower axis. 



3 One-dimensional toy model 



We can gain much understanding of the main features of the radio pulses in 
the time-domain with the aid of a very simple model. We show in this section 
that many of the characteristics of the pulse such as start-time, peak value of 
the electric field and duration in time, are mainly determined by the geometry 
of the system formed by the shower and the observer. In the model we assume 
a vertical atmospheric shower in which particles propagate along the shower 
axis at the speed of light c, and we ignore their lateral spread. For the moment 
we assume a constant index of refraction n. We do not make any assumption 
on the radio emission mechanism, and we just account for the retarded time, 
i.e. the radio signal is emitted at a time t' and reaches the observer located 
at a distance r to the shower core at a later time t = t' + Atp, where Atp is 
given by the electromagnetic wave travel time between the emission point and 
the observer. The model can also be applied to understand the radio emission 
properties in dense media 



In Fig. [2] we show a sketch of the model. The origin of time t' = is arbitrarily 
fixed at the time at which the shower front reaches ground. The emission time 
at a height h above the ground is t' = —h/c, the propagation time to an 
antenna on the ground at a distance r from the shower core is nR/c, where 
R = y/h"^ + r^, and hence the arrival time of the radio emission at the antenna 
is: 



' = 1 <»> 

The ID model although simple is useful to understand at least in a qualitative 
way the prominent role played by the geometry in the behavior of the electric 
pulse signal. 



3.1 Start-time of the electric field pulse 



The start-time time of the electric field pulse is given by the minimum value 
of t. By doing dt/dh = 0, one can find the height of the shower seen first by 
the observer: 



''start 



y^?^! 



(10) 



Substituting hstart in Eq- (E]), we obtain the time at which the observer sees 
the onset of the radio pulse: 



tstart = -Vn2 - 1 (11) 

C 

For n > 1, tstart is linear in r and the observer sees an increasing delay in the 
start time of the pulse as the distance to the shower core increases. 

In Fig. [3] we plot the relation given in Eq. (Q between the observer time t 
and the time at source t' or equivalently the shower depth (see also [2D]|). The 
relation is shown at different distances r to the shower core and assuming 
different values of n including the more realistic n varying with height (see 
section [273ll . One can see the increase of tstart with r when n > 1 as predicted 
by Eq. ( iTTl) . Also, the angle 6 between the observer and the shower axis (see 
Fig. [2]) is given by tan 6' = r/h. When h = hstart as given in Eq. (ITOl) it is 
straightforward to show that: 



tanOstart = y/n"^ - 1 = ta.n6c (12) 

and the observer sees first a shower height hgtart at time tstart with an angle 
equal to the Cherenkov angle. A similar analysis with similar results to the 
ones described above was developed independently and in parallel with ours 
by de Vries, Scholten and Werner [H], and in dense media in 



Due to relativistic effects associated to the speed of the shower (assumed to 
be c) being larger than the speed of the radio waves {c/n), at a fixed detection 
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Fig. 3. Relation in Eq. ([9]) between the observer time t and the time at the source t' 
or equivalently as depicted in the plot the depth of shower development. The relation 
is valid for vertical showers and is shown for observers at different distances r to the 
shower core and for different values of the index of refraction n. Also shown is the 
longitudinal profile of a 100 PeV proton shower (from a Gaisser-Hillas function) in 
arbitrary units, which gives an idea of the relative number of particles at different 
depths. 
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time td an observer at a distance r from the shower core may see two different 
stages of shower development (two different heights h±) simultaneously. This 
can be seen in Fig. 121 Solving the quadratic Eq. i^ for h with td constant we 
obtain: 



ctrf±ncJt^-r2(n2-l)/c2 
^^ = W^) ^ ^ 

Two different real solutions exist if ra > 1 and the argument of the square 
root is positive, which is equivalent to the condition td > t start- This apparent 
violation of causality is simply a relativistic effect due to the particles traveling 
faster that the speed of the pulse in the medium. When n = 1 the shower is 
seen by any antenna away from the core in a "causal" way, from beginning 
to end, as can be seen in Fig. [31 provided the observer is not located on the 
shower axis (when n = 1 and r = the whole shower is seen at the same 
instant of time). 

When r is sufficiently large or n = 1 (see Eq. (TTOj) ). hstart can become larger 
than the actual height at which the shower starts developing, namely the 
height /iq where the first interaction occurs. The emission from heights above 
the first interaction point is due to the primary particle only and can be 
neglected. In that case, the observer would start to see the onset of the pulse 
at a time corresponding to the height ho: 



— (14) 



The height /iq can be considered a physical limit on hstart ■, which corresponds 
to the depth Xq of the first interaction. The corresponding start times to can 
be read from Fig. [3] by truncating the curves in Fig. [31 at the minimum depth 
Xq. So, for any given ra, there is a critical distance Vcrit, at which hstart = ho: 



crit 



hoVn^ - 1 (15) 



For n > 1, as we move away from the shower axis the start time first increases 
linearly with r (following Eq. ( ITTj) ) until r = Vcrit- For larger r it follows the 
non-linear behavior of Eq. ([H]). For n = 1, the effective t start is always given 
by to in Eq. (JH]), regardless of the distance r to the core. 

Regarding the relativistic effects described above, when n = 1, r^rit = and 
the shower is seen in a causal way by any observer, as stated before. For 
n > 1, observers at distances greater than Vcrit will also see the shower in a 
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"causal" way, starting at Hq , and in this model they never see the shower at 
the Cherenkov angle. 



3.2 Peak value of electric field pulse 



The peak value of the radio pulse is mainly determined by three factors, 
namely, the distance from the shower to the observer, the number of charged 
particles in the shower, and also and in a very important way by geometrical 
effects associated to a time compression factor. The relation between the ob- 
server time and the shower height, given by Eq. i^ and shown in Fig. [31 is non 
linear. One can consider a given observer time interval and obtain from this 
relation the interval of the shower development that contributes to it. Clearly, 
the portion of shower development that contributes will be large when the 
slope of this relation is small. We can define a compression factor /^ taking 
the derivative of Eq. (^ with respect to h: 



fc 



dt 



dh 




[ns m '^] (16) 



The absolute value accounts for the fact that the derivative of t with respect 
to h changes sign when h = hgtart (i-e. t = t start) i corresponding to a reversal 
in the time sequence of the shower as seen by the observer. The inverse of 
the compression (or Doppler) factor /c can be thought of as a measure of 
the compression in time [32]. Qualitatively, a small value of fc implies that 
the emission from a relatively large portion of the shower contributes to the 
pulse in a relatively small interval of observer's time 5t, increasing the pulse 
with respect to other cases in which the factor /^ is larger. Moreover, when 
h — !■ hstart (i-e. t — )■ tstart) then /c — !■ and 5t -^ 0. This is related to 
the fact that the observer sees the shower at t = tgtart with an angle equal 
to the Cherenkov angle, as shown in Eq. (TT2|) . In fact, since cos 6' = h/R 
it is straightforward to show that the factor in Eq. (TT6l) is proportional to 
1(1 — 72/3 cos 6*) I . As a consequence of fc — )■ 0, the factor (1 — nf3 cos 9) in the 
denominator of Eq. ([3]) goes to zero, tending to enhance the peak value of the 
pulse, seen at the Cherenkov anglqlJ 

It is important to note that although a small fc induces a very big effect in 



'^ This singularity is not a problem for our numerical calculations for several reasons: 
For the singularity to actually enter the numerical calculation, fc would have to 
vanish (within the precision of our code) at precisely the middle of the track (see 
Fig. [T]). This (almost) never happens, but if it does, the code uses the limit of Eq. ([3]) 
for /c — 7" [25]. Also note that just like in reality (and Monte Carlo too) there is 
a finite time resolution (we use 6t = 0.5 ns in this work), which spreads very high 
(and short) vector potential peaks over the whole time bin. 
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the pulse in the ID model where all particles follow the shower axis, in a real 
shower this effect will be less important due to the lateral and angular spread 
of the particles in the shower. The lateral and angular spread will change the 
observer angle and increase the compression factor at hstart with respect to 
what the ID model predicts (/c ~ 0), since fc oc — cos^. Furthermore, the 
lateral spread of the particles in a real shower will also change the distance to 
the observer, smearing the arrival time of the signal. 

As stated before, the peak value of the pulse is also determined by other factors 
besides fc- There is a non-trivial interplay between the length and height of 
the part of the shower seen with a small compression factor, the number of 
particles in that region of the shower development, and the distance R to the 
observer. A clear illustration of this interplay can be seen in Fig. HI where 
we plot fc as a function of h along with a Gaisser-Hillas parameterization of 
shower development for a 10^^ eV shower. When n = 1.0003 (left panel) the 
largest compression in time applies to the region around shower maximum 
for an observer at r = 100 m, while for an observer at r = 50 m the largest 
compression applies to a portion of the shower with fewer particles, below 
shower maximum. 

One could think that the peak of the pulse will always drop as r increases and 
the observer gets further away from shower axis, but the factor fc can slow 
down this trend (e.g. between r = 50 m and r = 100 m in Fig.H]) or even reverse 
it. Moreover, as r increases, the width of the peak in fc becomes larger (this 
can be appreciated by inspection of the widths at values of fc = 10~^ ns/m in 
Fig. H]) and this implies that a larger portion of the shower contributes quasi- 
simultaneously to the observed emission. On the other hand, as r increases so 
does the distance from the emission point to the antenna, tending to decrease 
the signal. One can see that even in this simple ID model the interplay of the 
various relevant variables is already very complicated and somewhat counter- 
intuitive. For observers at larger distances to the core (e.g. r = 400 m in the 
bottom panel of Fig. [3]), the largest compression is achieved only at very high 
altitudes where the number of particles in the shower is much smaller than at 
shower maximum, and thus the peak value of the pulse is much smaller than 
that seen by an observer at r = 50 m (bottom panel of Fig. [3] and Fig. H]). 



3.3 Time duration of the electric field pulse 



The ID model also allows to understand the time duration of the electric pulse. 
The shower arrives at ground at t' = and the signal from the shower when 
it reaches the ground arrives at the observer at a time tg = nr/c. Assuming 
that the emission is only due to the shower front, the time duration At of the 
pulse, as seen by the observer, is then given by: 
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Fig. 4. Left Panel: Compression factor fc as defined in Eq. (I16p for n = 1.0003 (i.e. 
the absolute value of the derivative of the curves labeled n = 1.0003 depicted in 
Fig. [3|) as a function of altitude h for observers at different distances r to the shower 
axis. Also shown is the number of charged particles in a 10^^ eV proton shower as 
given by a Gaisser-Hillas (GH) parameterization of the longitudinal shower devel- 
opment (right axis). Right Panel: Same as left, but for n = n{h). 



At 



T^g T^start 



- [n — yri? 



(17) 



where t start is given by Eq. f lTT|) . This equation is only valid for r greater 
than a few meterqiJ and r < Tcrit (see Eq JTSj) . If r > Tcrit-, t start should be 
replaced by to, as defined in Eq. ^\M . but this does not change significantly 
the behavior of At, which stays approximately linear with r. So the observers 
see an increasingly wider pulse in time the farther they are from the shower 
core. Another estimate for the variation of the pulse width with distance can 
be found in [331. 



This approach is an oversimplification as it assumes that the shower still con- 
tributes to the pulse as it reaches ground level, so it needs to be modified for 
very inclined showers. In addition it also neglects the delays of the particles 
that lag behind the shower front and also contribute to the pulse, making it 
wider. Still the simple expressions obtained can be quite useful to interpret 
the results from the full simulation, shown in the next section. 



3.4 Dependence on refractive index 



The ID model also allows to study the dependence of the pulse properties on 
the index of refraction. For the purpose of understanding the relevance of n 
on the pulse, we first assume n to be constant with height. Using Eq. (1111) 
it is straightforward to deduce that the start-time of the pulse increases with 
n because the propagation time of the signal from the height of emission to 
the observer increases with n. Furthermore, the time duration of the pulse 



For very small distances r < ho{n — l)/2n, then to > tg and At = to — t 



start' 
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also depends on n. As can be seen in Eq. (fT7|) . the time duration of the pulse 
decreases slightly with increasing n. 

Another dependence of the pulse properties on n arises because the factor /^ 
is also dependent on n. An observer at a fixed given distance r from the core 
will see different parts of the same shower with a small fc (i.e. with a large 
compression in time) for different values of n, since hstart (which defines the 
altitude at which fc = and the shower is seen with 6 = 9c) decreases with 
n. To illustrate this we show in Fig. O the factor /^ as a function of h for an 
observer at r = 100 m for several values of n. One can see that as n decreases 
the height at which fc = (i.e. hstart) increases, leading to different parts of 
the shower being largely compressed in time. The effect on pulse height will 
depend on the number of particles at hstart^ which in turn depends on the 
distance r to the observer. For r < 100 m (not shown in Fig. [S]), the highly 
compressed part of the shower for n = 1.0003 is below shower maximum. As 
n decreases, the compressed part of the shower moves towards the maximum, 
tending to increase the pulse. On the other hand, for r > 100 m, the height 
hstart is above the maximum, and a decrease in n will move the compressed 
part of the shower further away from the maximum, tending to decrease the 
pulse height. For large values of r, hstart is in a region with very few particles, 
decreasing the effect of a low fc in pulse height. Furthermore, in this part of 
the shower the number of particles increases only very slowly with decreasing 
h, since the density at these altitudes is low, and the effects of changes in n 
on pulse height become less important. 

As discussed before, relativistic effects such as seeing the later parts of the 
shower before the earlier ones and seeing two parts of the shower simultane- 
ously can only be observed at distances r < rcrit from the core. Since Vcru 
increases with n, observers further from the core will also see these effects as 
n increases. 



4 ZHAireS simulations 



In this section we present the electric field as obtained in ZHAireS simulations 
of atmospheric showers. In particular in this section we will obtain the behav- 
ior of different features of the time pulse with r and refractive index, but this 
time in realistic simulations that take into account the full complexity of atmo- 
spheric showers. We will see that the behavior follows the qualitative behavior 
obtained with the ID model. At the end of this section we also show some 
results of the electric field calculated in the frequency domain and compare 
them with Fourier transforms of the time pulses. 
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Fig. 5. Same as Fig. [U but for an observer at r = 100 m and several values of n, 
including a model for the variation of n with altitude n{h) (see Section [2. 3|) . 

4-1 Dependence of the electric field pulse on distance to shower axis 



In Fig. |6] we show the East-West (EW) electric field component as a function 
of time obtained from simulations of 10^^ eV vertical atmospheric showers 
induced by protons. We used a horizontal magnetic field \B\ = 23 fiT pointing 
north for this particular simulation. The electric field is shown at different 
distances r northwards from the shower core. As predicted by the ID model 
in Eqs. flTTl) and flT7|) . it can be seen that the start-time of the field and the 
duration in time (at least the easily visible positive part of the pulse) both 
increase with r for the observers shown in the plot. 

The peak of the electric field decreases with r for an observer North of the 
shower core, but not in a linear manner. In fact it decreases roughly as r from 
50 m to 100 m, but faster than r from 100 m to 150 m, as can be seen in Fig. O 
As explained above, there is a non-trivial interplay between the distance from 
the shower to the observer, the factor fc and the number of particles in the 
region of the shower seen with a large compression in time, illustrated in the 
right panel of Fig. HI The number of particles at hgtart increases only slightly 
from r = 50 m to r = 100 m, but decreases by an order of magnitude between 
50 m and 150 m. Since the distance R from emission point to observer increases 
roughly as i? ~ r, it is twice (three times) as big at 100 m (150 m) than at 
50 m. This explains why the height of the peak at 100 m is about half the 
height at 50 m, while it is much smaller at 150 m. We have also observed in 
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Fig. 6. EW component of the electric field at r = 50, 100, 150 and 400 m northwards 
of the shower core as obtained in ZHAireS simulations of a 10^^ eV proton-induced 
vertical atmospheric shower. These simulations were performed with the exponential 
model for the variation with altitude of the refractive index. Note that the signal 
at r = 400 m is arbitrarily multiplied by a factor 20 for plotting purposes. 

our simulations that the relative height of the pulse at different distances also 
depends on the direction of the observer w.r.t. the shower core. In Fig. [7] we 
show the emission from a similar shower as in Fig. [6l but for observers East of 
the core. For an observer East of the shower core, the peak amplitude appears 
to be maximal at around r = 100 m, a result in principle compatible with |?T] . 
Also, the decrease in pulse height with r seems to be slower in Fig. [71 when 
compared to the observer North of the core (Fig. [6]). A similar dependence on 
antenna position can also be seen for the spectra at 1 MHz, shown in FiglTSl 
This difference is, in part, due to the interference of the geomagnetic and 
Askaryan components of the emission, as will be discussed in section [5l 



4-2 Dependence of the electric field pulse on the refractive index 



In Fig. [H] we show the EW component of the electric field as a function of time 
obtained from simulations of 10^^ eV vertical atmospheric showers induced 
by protons, using a magnetic field \B\ = 23 /iT with an inclination of —37° 
and a declination of 0°. The calculations were performed with two constant 
refractive indices, namely n = 1.0 and n = 1.0003, as well as with a refractive 
index varying with altitude according to Eq. ([7]). 
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Fig. 7. Same as Fig. [6l but for observers east of the core. 



Firstly, the larger the refractive index the later the pulse starts as obtained 
with the ID model above (see Eq. [TT1) . The larger the n the smaller the prop- 
agation speed, and so the signal reaches the observer at later times. For the 
refractive index changing with altitude, which varies between n = 1.0 and 
n = 1.0003, the start-time falls between the two obtained for the two constant 
refractive indices, as expected. 



For observers close to the shower axis (r = 100 m in the left panel of Fig. |8]), 
small changes in n are responsible for large changes in peak height and width. 
This can be qualitatively understood, as described in Sections 13.21 and 13. 4[ in 
terms of the non-trivial interplay between the various geometrical factors. In 
Fig. m and E] one can see that for r = 100 m the number of particles in the low 
fc region changes only slightly from ~ 70 ■ 10^ when n = 1.0003 to ~ 60 ■ 10^ 
when n = n{h), while the length of the shower seen with fc < 10~^ doubles. For 
observers further away from the core (r = 400 m in the right panel of Fig. |8]), 
the effect of the varying n in pulse height is much less pronounced, since /c ~ 
only above the shower, as can be seen in the bottom panel of Fig. 121 This will 
make the compression factor in the shower region very similar for different 
values of n, as illustrated by the fact that the curve labeled r = 400 m in 
Fig. m changes only very slightly from the left panel (n = 1.0003) to the right 
panel (n = n{h)). 
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Fig. 8. EW component of the electric field at r = 100 m (left panel) and r = 400 m 
(right panel) eastwards from the core for a vertical 10^^ eV proton-induced shower as 
obtained in ZHAireS simulations. The calculation was done using several refractive 
indexes. 

4-3 Dependence of the electric field pulse on the shower zenithal angle 



The features of the electric field at a fixed distance to the shower axis are 
expected to depend strongly on the shower zenithal angle. The main reason 
for this is that the shape of the curve relating the observer time and the source 
depth - which largely determines the characteristics of the pulses as explained 
above - depends strongly on the geometry of the system formed by the observer 
and the shower. The electric field from a certain region of the shower depends 
on the angle 6i (with which the observer sees that region) through the factor 
|1 — nP cos 6i\ in Eq. ([3]), on the number of particles Ni in the region, and 
on its distance to the observer, Ri. For 6i close to the Cherenkov angle, the 
emission is enlarged with respect to other angles, since \1 — n cos di \ is small. 

The interplay between |1 — ncosdi\ (or equivalently the factor /c) and N^ is 
shown in Fig. [9] for a vertical (zenithal angle 6 = 0°) and an inclined shower 
{6 = 50°). The curves in Fig. [9] were obtained for an observer at a distance 
r = 400 m (in the early part of the inclined shower) and n = n{h) with a 
modified ID model, similar to the one developed in Section [3], but which can 
handle inclined showers. Fig. [9] suggests that inclined showers produce larger 
signals than vertical showers for distances greater than a couple of hundred 
meters from the core. The main reason for this is that at larger distances to the 
shower axis, inclined showers are viewed with angles closer to the Cherenkov 
angle than vertical showers, and so the angular factor 1 1 — n(3 cos 6i \ is smaller, 
boosting the emission and the net signal from inclined showers. A similar 
conclusion was also reached in 



Using ZHAireS, we investigated the dependence of the field on zenithal angle 
as a function of distance to shower axis, in a full Monte Carlo simulation of 
the shower development in air with a realistic model for the refractive index. 
We compared the signals obtained from 10^^ eV proton-induced showers with 



19 



|1-ncosa|-e=50' 



|i-ncosej -e=o° 

GH - e=50° 



80 



0) 

u 
50 ■■? 




2000 4000 6000 8000 10000 1200014000 1600018000 

h(m) 



Fig. 9. Geometrical angular term |1— ncos^il as a function of altitude /i for a vertical 
shower (dashed blue line) and a shower with zenithal angle 6 = 50° (solid red line). 
Also shown are Gaisser-Hillas parameterizations of the number of particles (linear 
right axis) as a function of height for a lOOPeV vertical shower (thick dashed blue 
line) and an equivalent one with 6 = 50° (thick solid red line). The calculations 
were done numerically for n = n[h) and r = 400 m (early part of inclined shower), 
using an ID model that allows non-vertical showers, similar to the one described 
for vertical showers. 

^ = 0° and 9 = 50°. We found that the results obtained with ZHAireS are 
in qualitative agreement with our ID model and with the calculations in [3^ . 
In Fig. [To] we compare the East- West components of the electric field for 
the vertical and inclined showers for an antenna at r = 100 m (left) and 
r = 800 m (right) located in the North-East of the shower core. Clearly at 
large distances to the core the peak value of the electric field in the inclined 
shower is much larger than that of the vertical shower, a factor ~ 5 in this 
particular simulation. This trend of an increased non-vertical shower signal at 
large distances from the core is present in all directions with respect to the 
shower core. Note also the shift in the start-time of the pulse in the inclined 
shower w.r.t. the vertical one, due to the different geometry of both showers. 

We have also found that the signal in non-vertical showers at large distances 
from the core is very sensitive to the decrease of the refractive index with 
altitude. In Fig. [TT] we plot the electric field for the same showers as in Fig. 
[T0| but instead of the model of variable refractive index with altitude, the fields 
shown were calculated using a constant refractive index n = 1.000325. It can 
be clearly seen in the left panel of Fig. [11] that for the antenna at r = 100 m 
North-East from the core, the use of a constant n decreased both the vertical 
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Fig. 10. Comparison between the electric field as a function of time as obtained in 
ZHAireS simulations of 10^^ eV proton-induced showers with zenithal angle 9 = 0° 
(blue solid) and 6 = 50° (red dashed), for antennas at r = 100 m (left) and r = 800 m 
(right) NE of the shower core on the ground. The inclined shower comes from 
the SE. The simulations were done using a variable refractive index model for the 
atmosphere. 
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Fig. 11. Same as Fig. [TOl but using a constant refractive index n 
whole atmosphere. 



1.00035 for the 



and non- vertical EW field component w.r.t. the variable n, while at r = 800 m 
(right panel), the signal in the vertical shower stayed practically unchanged 
and the peak in the non- vertical shower more than doubled in value. Again this 
behavior can be traced back to the geometrical dependence of the interplay 
between the various key elements in the calculation of the electric field. 



4-4 Radio pulses in the frequency domain 



In Fig. [12] we compare the results of ZHAireS frequency domain calculations 
with Fast Fourier Transforms (FFT) of the calculated time domain signals for 
the same shower. The convention used for the normalization of the Fourier 
transform is in Eq. ([5]). One can see that the agreement between the spectra 
and the FFT is very good up to high frequencies, at which the width of the 
positive peak in the time domain becomes important, since it is the smallest 
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Fig. 12. Comparison between the emission spectrum calculated by ZHAireS in the 
frequency domain, and the FFT of the time domain signal for the same 100 PeV 
vertical shower. 

large scale structure of the time pulse. More quantitatively, above v ~ 1/AT, 
where AT is the characteristic width of the peak of the pulse in time, the 
FFT starts fluctuating. For m (400 m), the peak width is around 5 ns 
(50 ns), leading to large fluctuations starting at around 200 MHz (20 MHz)[£j. 
This high frequency part of the spectrum is largely incoherent and sensitive to 
shower to shower fluctuations and to the thinning level used in the simulation. 
These fluctuations are partly physical (incoherence level, as will be discussed 
below) and partly unphysical (bin size effects, thinning, FFT numerical er- 
ror, etc). Also, in the direct frequency domain calculations they are slightly 
smaller than in the Fourier-transform of the time pulsa ^ I. This means that to 
study the spectrum at high frequencies it is best to calculate shower emission 
directly in the frequency domain, reducing the unphysical contributions to the 
fluctuations. 



The spectrum of the radio emission depends on the distance r from the core 
on the ground to the antenna. As can be seen in figure [T21 the frequency at 
which the Fourier components reach a maximum decreases from ~ 30 MHz 
at r = m to ~ 3 MHz at 400 m. This cutoff frequency is related to the 
time duration of the pulse for a specific observer and serves as a boundary 



^ Similar fluctuations can also be seen in [17] • 

^ Note that due to computing time issues, the frequency bin used in the frequency 
domain simulations shown in Fig. [12] is much wider than the frequency bin used in 
the FFT, making the fluctuations appear much smaller than they really are. 
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between the fully coherent and incoherent parts of the spectrum. Above this 
frequency, the emission is subject to shower to shower fluctuations. At even 
higher frequencies, effects due to thinning in the simulation become important, 
since the spectrum becomes sensitive to the fine details of the shower, even 
down to single particles for the highest frequencies (not shown). A 2D model 
that explains this behavior in the Fraunhofer regime, but which is still relevant, 
can be found in [53]. 

The fact that the cutoff frequency of the emission decreases with distance to 
the core has implications on the measured lateral distribution of the signal. 
Depending on the measured frequency range, the Fourier components of the 
signal will drop at different rates as the distance increases. In Fig. [13] we show 
the average total field as a function of distance to the core at 1 MHz (left) 
and at 60 MHz (right), along with their RMS, calculated from simulations 
of 10 proton showers of 100 PeV. At 60 MHz, i.e. in the incoherent part of 
the spectrum, the field decreases much faster with increasing distance to the 
core compared to the coherent part of the spectrum at 1 MHz. This behavior 
can be relevant for the design of radio detector arrays. An explanation of the 
reason why the largest (smallest) pulses are received at antennas located in 
the East (West) of the shower is given in the next section. 

At the fully coherent 1 MHz frequency (Fig. [T3]left), one can also see that a 
maximum appears at around r ~ 100 m. A similar maximum was reported 
in [H], but for full bandwidth pulse height. Furthermore, Fig. [13] suggests 
that this maximum is dependant not only on frequency, but also on observer 
direction w.r.t. the shower core, disappearing for antennas to the west. Our full 
bandwidth pulses (Figs. [H] and [7|) show a similar dependence of the maximum 
of the emission with observer direction. This suggests that the mechanism 
responsible for this maximum is weaker than the interference effect between 
the geomagnetic and Askaryan components of the emission, which is discussed 
in the next section. 



5 Polarization properties of the signal and asymmetries 



The polarization of the electric field depends on the relative importance of the 
two main mechanisms thought to be responsible for the radio emission [20]. 
Cherenkov radio emission due to the excess of electrons in the shower (Askaryan 
effect) has a characteristic axial polarization in relation to the shower axis. On 
the other hand, the polarization of the emission due to geomagnetic follows 
—V X B, and depends on the direction v of the charged particles and the geo- 
magnetic field B. In contrast to the charge excess mechanism, the geomagnetic 
emission of electrons and positrons add up since they have opposite charges 
but their trajectories curve in opposite directions (see Section [TT]) . In the top 
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Fig. 13. Average electric field vs distance to core for 10 vertical showers of 100 PeV 
simulated with ZHAireS at 1 MHz (left) and 60 MHz (right). The observers are 
placed at positions N, S, E and W of the shower core. Note that for observers N 
and S the fields are similar and lie on top of each other in the figure. Also shown as 
points are the RMS of the electric field. 

panel of Fig.[Tl]we show schematically the projection of the polarization on the 
ground induced by each of the emission mechanisms for a vertical shower [17] . 
A horizontal magnetic field pointing north (N) is assumed. 

The interplay between the different polarizations makes the projection of the 
net field on ground asymmetric with respect to the shower core [l5] . This can 
be easily understood from the sketch at the top of Fig. [TH Since the Askaryan 
and geomagnetic polarizations point in the same direction for observers East 
of the core, and in opposite directions for observers to the west, we expect the 
EW component of the electric field (Eew) to be larger in an antenna eastwards 
of the shower core than in an antenna West of the core, at the same distance. 
We also expect Eew to be of the same order northwards and southwards from 
the shower core. Another prediction is that the NS component of the electric 
field (Ens) should be close to zero for observers along the East- West direction 
and largest for observers along the North-South one ^7j. As discussed in |15] . 
as the distance r from the antenna to the shower axis increases, the relative 
contribution to the electric field of the geomagnetic mechanism decreases and 
as a consequence the Askaryan mechanism accounts for a larger fraction of 
the emission. The pattern of the electric field at ground level is then expected 
to recover a symmetric behavior with respect to the shower core for very large 
distances, when the Askaryan contribution dominates. 

To investigate these expectations we have obtained the Eew and Ejys com- 
ponents of the field, at a frequency of 60 MHz, from ZHAireS simulations of 
vertical proton showers with energy 10^'^ eV. These components are shown in 
Fig. [TH In the left panel we show E^s as a function of the position on the 
ground, while in the right panel we show the corresponding Eew- One can 
see the expected EW asymmetry in the Eew component, with larger fields 
to the east of the core, while it is approximately the same North and South 
of the shower core. The E^s component (left panel) is mainly due to the 
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Fig. 14. (Color online) Top: Sketch of the projection on the ground of the electric 
field induced by the Askaryan and geomagnetic emission mechanisms in a vertical 
shower 9 = 0° (see also [17])- Bottom: Left panel: North-South component of the 
electric field E^s as a function of the position around the shower core as obtained 
in ZHAireS simulations of vertical showers induced by protons of energy W^'^eV. 
Right panel: East- West component Eew of the electric field obtained in the same 
simulations. The color scale indicates the magnitude of the components of the field 
in (V/m/MHz), note the different scale in the left and right panels. 



Askaryan mechanism, and hence it is largest along the NS direction, while it 
gradually decreases as the observer moves to the EW direction because nei- 
ther the Askaryan nor the geomagnetic mechanism induce a significant Ens 
component along observers in the EW direction (see sketch in Fig. E 



Comparing Ens and Eew for observers along the NS axis, one can clearly 
see that at distances relatively close to the shower core (r < 150 m), the 
NS component of the field is a factor ~ 4 smaller than the EW component, 
because the geomagnetic mechanism dominates the emission close to the core, 
while this factor tends to diminish at larger distances (see also Fig. [T3ll . 
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5.1 B = vs B ^ 0: Separation of the geomagnetic and Askaryan compo- 
nents. 



In order to disentangle the geomagnetic and Askaryan components of the ra- 
dio emission, we also simulated showers turning off the geomagnetic field. As 
discussed in section [5l there is an EW asymmetry in the signal strength due to 
an interference effect between the polarizations of the different emission mech- 
anisms, making the signal larger to the East. In fig. [15] we compare the EW 
components (positive to the East) at 100 m (top) and 400 m (bottom) West 
(left) and East (right) of the core, obtained with simulations with and with- 
out the geomagnetic field. The simulation with the magnetic field on shows 
the net field due to both emission mechanisms, while the simulations without 
the magnetic field has only the Askaryan component. One can see that, as 
expected, the polarization of the Askaryan component changes direction be- 
tween antennas East and West of the core (dashed red lines of fig. [15]), and 
this causes the difference in the peak height of the net field (solid blue lines). 
Also, the height of the pure Askaryan peak (dashed red) is roughly half the 
difference between the E and W peaks of the net field (solid blue), as ex- 
pected [T7] (see also the sketch in Fig. [H]) . At larger distances (e.g. at 400 m, 
at the bottom of fig JT^ , the ratio between the Askaryan and geomagnetic 
components gets higher. At larger distances, the Askaryan mechanism starts 
to dominate the emission. 

To further investigate the polarization of the emission we use a parameter TZp, 
as defined in [H], which is sensitive to the polarization of the electric field. In 
the case of a horizontal (parallel to ground) magnetic field pointing north, the 
polarization vector of the geomagnetic contribution to the electric field points 
East, 



and TZp is given by: 




7^p = 


Y.t Eew ■ Ens 
Et {E%w + E%s) 



(18) 

where the sum runs over all bins in time having a non-zero electric field at the 
antenna. 

From the sketch in the top of Fig. [H] we expect TZp to vary as a function of the 
azimuthal angle on the ground, defined so that = for an antenna in the 
East direction and (p = 90° for an antenna at the North. TZp = when either 
Eew ot Ejys equal zero. We then expect that if the polarization were only 
due to the geomagnetic contribution (the unphysical case of no charge excess 
in the shower), the NS polarization would always be zero and so would TZp. 
If at a certain distance r to the shower core the Askaryan and geomagnetic 
contributions are important, then TZp ~ at = 0° and = 180°, because 
Ens ~ along those two directions [S]. However, in contrast to what was 
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Fig. 15. Comparison between the EW component of the field calculated with and 
without the geomagnetic field at a 100 m (top) and 400 m (bottom) West (left) and 
East (right) of the core. 



expected in [44], if the geomagnetic component is absent then we clearly expect 
7^p ~ at = 0°, 90°, 180° and 270°, because either Eew ~ or Ens ~ 
along those directions. This means that in the absence of a magnetic field, 
TZp should exhibit a 180° periodicity in azimuthal angle (p. We checked this 
in our ZHAireS simulations by switching off the magnetic field. The result is 
plotted in Fig. [16] where one can see that TZp exhibits a periodicity of 180° in 
0, regardless of the distance to the observer. 



In contrast, in Fig. [17] we show TZp as a function of as obtained in ZHAireS 
simulations of 10 vertical proton showers of energy 10^^ eV, but in this case 
with the geomagnetic field on. One can see that for distances closer to the core 
(r < 300 m), the TZp parameter reaches zero only at = 0° and = 180°, as 
expected when both the geomagnetic and Askaryan components are impor- 
tant, while at larger distances from the core, the period of TZp vs changes 
from 360° to ~ 180°. This further confirms that at large distances to the core 
(r > 600 m), the Askaryan mechanism dominates the full bandwidth emis- 
sion. A similar behavior was reported in |13] for REAS3 simulations convoluted 
with the detector response, but this periodicity change in TZp was interpreted 
as a signature of a dipole polarized field instead of a dominant charge excess 
contribution. In the same paper |44j, MGMR simulations did not exhibit this 
behavior. 
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Fig. 16. ZHAireS code predictions for the polarization sensitive parameter TZp (as 
defined in Eq. (jlSp ). as a function of the azimuthal angle (p (see top of Fig. [T^ for 
10 vertical proton showers of energy 10^^ eV. The magnetic field was switched off 
in the simulations. The TZp parameter is plotted for antennas close to the shower 
core (blue circles: r < 300 m), and farther from it (red stars: r > 600 m). 

If we assume a perfect shower symmetry in 0, the modulus E'^^ of the electric 
field due to Askaryan mechanism would be the same for antennas at the same 
distance from the core, and the EW and NS components could be written as: 



E^^ = -E^'^ cos (f), E'^). = -E^^ sin , 



(19) 



If we further assume no time dependence of the polarization, and take the 
numerator of Eq. (fTSll : 

7^J'^oc(E^^)%in20 , (20) 

which leads to TZp^ with a period of 180° in 0, as shown in Fig. [161 

The field with both the Askaryan E'~"^ and geomagnetic E^'^° contributions 
can be written as (see top of Fig. [1] 



Eew = - {E^"° + E^^ cos 0) , E_ 



leading to: 



NS 



Ch\2, 



E^^sin. 



Up oc E^"" E^^sin0 + (E^^) sin 20/2 



(21) 
(22) 
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Fig. 17. ZHAireS code predictions for the polarization sensitive parameter TZp (as 
defined in Eq. (jlSp ). as a function of the azimuthal angle (p for 10 vertical proton 
showers of energy 10^^ eV, with a horizontal magnetic field pointing north. The TZp 
parameter is plotted for antennas close to the shower core (blue circles: r < 300 m) 
and farther from it (red stars: r > 600 m). 

The equation above predicts that when the Askaryan and geomagnetic mecha- 
nisms compete, TZp is proportional to a rather complicated function of (p whose 
shape depends on the relative values of E*^^ and E^^°. 

In Fig. [18] we show the numerator of TZp, given by the product Eew ■ E^s, as 
a function of obtained in ZHAireS simulations of a vertical proton shower 
of 10^^ eV, for antennas at r = 200 m (left) and r = 1200 m (right) from the 
shower core. We then fitted Eq. (122!) to these simulations with only E'^^ and 
E^^° as free parameters. The fit based on the simple model above reproduces 
remarkably well the behavior with obtained in the simulations, with a ratio 
^geo/^ch _ ^74 at r = 200 m, and E^^°/E^^ ~ 0.51 at r = 1200 m, de- 
creasing with r as the Askaryan component is expected to dominate at large 
distances to the core. It is interesting to note that in the ZHAireS simulation 
for antennas south of the core, the ratio between the peaks of Eew (expected 
to be purely geomagnetic) and E^s (expected to be purely Askaryan) de- 
creases from 2.32 at r = 200 m to 0.57 at r = 1200 m, for this same shower. It 
is also interesting to see in Fig. [18] how the parameter TZp gradually changes 
from having a 360° periodicity in to a 180° periodicity as the distance to the 
shower axis increases. We believe that similar, but more refined analysis meth- 
ods could be derived to separate the contributions of the different emission 
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Fig. 18. Points: The numerator of TZp in Eq. (ll8p as a function of (p for a vertical 
proton shower of 10^^ eV simulated with ZHAireS, for antennas at r = 200 m (left) 
and r = 1200 m (right) from the core. Solid lines: Fit of Eq. (j22p (obtained using a 
simple model of the polarization - see text for more details) to the simulated electric 
fields. 




(])+180 



Fig. 19. Sketch of the electric field in non vertical showers and its projection onto 
the ground plane. 

mechanisms to the electric field. 



5.2 Polarization of the electric field in non-vertical showers 



In the case of non-vertical showers, the projection on the ground of the po- 
larization of both, the geomagnetic and Askaryan components will depend on 
the azimuthal angle of the antenna position on the ground in relation to 
the core, as can be seen in schematic form in Fig. [191 Furthermore, the di- 
rection of the shower axis with respect to the geomagnetic field will change 
the direction of the Lorentz force acting on the charged particles. In vertical 
showers the force is always almost parallel to the ground, but in non-vertical 
showers the Lorentz force will only be horizontal if the plane defined by B and 
the shower axis is perpendicular to the ground (e.g. Fig. [20] top left). In most 
geometries there will be a vertical component to the Lorentz force (e.g. Fig. 
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top right). This will cause a dependence of the polarization and the signal 
asymmetries on the azimuthal angle of the shower, as can be seen in Fig. [201 
which shows the results of full ZHAireS simulations of the N-S (top), E-W 
(middle) and Z (bottom) components of the electric field at 60 MHz for 100 
PeV proton shower with 9 = 45° coming from the north (left) and the west 
(right), corresponding to the geometries shown schematically on the top left 
and top right of Fig. |20l respectively. Note that different scales were used for 
the field in the different panels. A horizontal magnetic field pointing north was 
used in the simulations. For the shower coming from the north (left), one can 
see that the main asymmetry is in the E-W component (left-middle), similar 
to the vertical shower case, while the N-S and Z components are smaller and 
very similar to each other, because the dominant geomagnetic contribution 
is horizontal. On the other hand, in the case of the shower coming from the 
West (right), there is a large asymmetry to the East on both, the 15- W and 
Z components, since in this particular geometry the dominant geomagnetic 
contribution makes an angle of ~ 45° with the horizontal, and thus should 
have very similar E-W and Z components. 

This dependence of the polarization on the azimuthal angle of the shower di- 
rection is relevant for studies trying to disentangle the contributions of the 
emission mechanisms using polarization. Since the inclination of the polariza- 
tion vector of the dominant geomagnetic emission with respect to the hori- 
zontal plane changes with the shower azimuthal angle, it may be important to 
also measure the vertical component of the net electric field, since for showers 
with 6 > 45°, the vertical component of the field can be even larger than the 
horizontal ones. 



6 Conclusions 



In this work we present predictions for the radio pulse emitted by extensive 
air showers. Our results are obtained using the ZHAireS Monte Carlo, an 
AIRES-based code that takes into account the full complexity of ultra-high 
energy cosmic-ray induced shower development in the atmosphere, and allows 
the calculation of the electric field in both the time and frequency domains 
based on the algorithms developed in [23112^ . Although our approach does 
not presuppose any a priori emission mechanism, our results confirm that 
the emission at radio frequencies can be understood as the superposition of 
radiation from the charge separation induced by the magnetic field of the 
Earth (geomagnetic effect), and that coming from the net excess negative 
charge evolving as the shower develops in the atmosphere (Askaryan effect). 

We have pointed out the relevance of the refractive index in the time structure 
and intensity of the radio pulses, especially at short distances to the shower 
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Fig. 20. Results of full ZHAireS simulations of the N-S (top), E-W (middle) and Z 
(bottom) components of the electric field at 60MHz for lOOPeV proton shower with 
6 = 45° coming from the north (left) and the west (right). See text for more details. 
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axis. Another interesting work [H], developed independently and in paral- 
lel with ours, deals with similar issues. The refractive index determines the 
angular distribution of the radiation at the emission point, as well as its prop- 
agation through the atmosphere and non trivial relativistic effects arise due to 
the refractive index n > 1. We have developed a simple 1-dimensional model 
to address the characteristics of the radio pulse, which qualitatively allows us 
to interpret the pulse height, width and its dependence on the refractive index 
and the distance from antenna to shower core j41j . The intensity of the radio 
pulses is typically highest for those observers that see a region of the shower 
containing a large number of particles (such as shower maximum) with an 
angle close to the Cherenkov angle. There is a non-trivial interplay between 
the distance from the emission point in the shower to the observation point, 
the angle between the particle direction and the observer, and the number of 
particles in the region of the emission, whose effects can only be accurately de- 
termined with Monte Carlo simulations such as those developed in this work. 
These key elements determine for instance that far from the shower axis in- 
clined showers typically induce a more intense radiation than vertical ones [32] 
as we show with our ZHAireS simulations. 

We have shown that the frequency at which the emission spectrum is max- 
imum, which is in the range ~ 3-30 MHz for < r < 400 m, decreases as 
the distance r from the antenna to the shower core increases (for r > 100 m 
see also [TT|17II46 ]). So most experiments, which typically are only sensitive 
to frequencies above 30 MHz, in fact only measure the incoherent part of the 
spectrum, even very close to the shower core. We also observed that the signal 
at higher frequencies decreases more rapidly with r than at the fully coher- 
ent 1 MHz range, where we observed a maximum in the emission at around 
r ~ 100 m, similar to the one reported in [H]. Our results also suggest that 
this maximum is dependant on frequency and on observer direction w.r.t. the 
shower core. We have also explored the polarization properties of the radia- 
tion, confirming the expectation that the radiation is mainly polarized in the 
opposite direction to the Lorentz force induced by the magnetic field of the 
Earth [28|29f30|l3T] . Far from the shower core, our simulations show that the 
polarization is compatible with the presence of a significant amount of radia- 
tion due to the Askaryan effect, as in [17]. We have also shown, using a very 
simple model, that as the Askaryan component of the emission becomes dom- 
inant over the geomagnetic one at larger distances to the core, the parameter 
TZp changes its periodicity on the azimuthal angle of the observer. When both, 
the geomagnetic and Askaryan components are significant, the periodicity is 
360°, but when the emission is dominated by the Askaryan effect, a 180° pe- 
riodicity can be clearly appreciated. In inclined showers we have stressed the 
role played by the relative orientation of the shower axis and the magnetic 
field on the polarization. A significant vertical component of the electric field 
arises in inclined showers which calls for detection instruments able to observe 
it. 
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